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momentum 



Abstract The zero-energy universal properties of scattering between a particle and a dimer that involves an 
identical particle are investigated for arbitrary scattering angular momenta. For this purpose, we derive an 
integral equation that generalises the Skorniakov - Ter-Martirosian equation to the case of non-zero angular 
momentum. As the mass ratio between the particles is varied, we find various scattering resonances that can 
be attributed to the appearance of universal trimers and Efimov trimers at the collisional threshold. 
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1 Introduction 

Few- and many-body systems at low energy with resonant short-range interactions are remarkable for their 
universal nature. In these systems, most of the details of the interaction potentials are unimportant except for a 
few parameters, which fully characterize the interaction between the particles. In many cases, these parameters 
reduce to a single one, the s-wave scattering length a, and all the physical properties are expressed in terms 
of the scattering length a. For example, the binding energy of two particles scales as a~ 2 , and the scattering 
properties between these dimers and other particles or other dimers can also be universally expressed in terms 
of a 01; SID]. When the s-wave scattering length becomes infinite, at some resonance point, one reaches the 
so-called unitarity regime, where the system has no intrinsic length scale and becomes completely universal 

a. 

In some other cases, however, extra parameters are needed. This is the case of three-body systems where 
the Efimov effect occurs ]5|; |6fl. The Efimov effect is characterized by the appearance of an infinite series 
of trimer bound states at unitarity. This effect results from the effective 3-body attraction between particles 
interacting through resonant short-range interactions. When this attraction is strong enough to overcome the 
repulsion that arises either from the centrifugal barrier due to rotation or from the Pauli exclusion principle 
of identical fermions, the particles can bind to form trimers called Efimov trimers. Because the particles 
are attracted to short distances, an extra short-range parameter is needed to completely determine their wave 
function. Thus, trimer binding energies and scattering properties in such systems depend on both the scattering 
length and the extra parameter. Efimov trimers are associated with remarkable features such as discrete scale 
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invariance of their spectrum, and have attracted a lot of interest since their recent experimental realizations 
with ultracold atoms §; S H M, O; E B • 

In general, the strength of the Efimov attraction depends on the mass ratios between particles. For example, 
in the case of two identical fermions and one different particle, which is resonantly interacting with the two 
fermions, the Efimov attraction is too weak for the Efimov effect to occur if the particles have the same 
masses. However, when the mass ratio between one of the two fermions and the third particle becomes larger 
than the critical value 13.6, and the system has one unit of total angular momentum, the Efimov attraction 
becomes dominant, and can support bound states, i.e. trimers formed of two fermions and one particle. 

Recently, Kartavtsev and Malykh found that trimers also exist for mass ratios smaller than the critical 
value [14]. Indeed, even though in this regime the centrifugal repulsion dominates at short distance, its com- 
bined effect with the weak Efimov attraction creates a shallow potential well that can support these trimers. 
Although these trimers are related to the Efimov physics, their nature differs from that of Efimov trimers. 
First, they do not exhibit the discrete scale invariance of the Efimov trimer spectrum, and appear only for 
positive scattering lengths. Second, they are universal states in the sense that they can be characterized only 
by the scattering length, unlike Efimov states which also depend on short-range physics. For the universal 
trimers, the centrifugal barrier at short distance prevents the particles from coming close to each other, so that 
the three-body recombination occurring at a short distance is also suppressed. Thus, these universal trimer 
states are much longer-lived than the Efimov trimers fl 1 511 - 

In this paper we investigate the general problem of the particle-dimer scattering in the zero-range limit, 
where two of the particles are identical (bosonic or fermionic) particles and a third particle is of a different 
kind. In particular, we look for the resonances in the particle-dimer scattering due to the formation of the 
universal trimers. For this purpose, we derive an integral equation to calculate the particle-dimer scattering 
lengths in various angular momentum channels. This equation is a generalization of the Skorniakov-Ter- 
Martirosian equation. With this integral equation, we calculate the particle-dimer scattering lengths for various 
angular momenta.We confirm that the scattering lengths diverge at the mass ratio where the universal trimers 
were predicted to appear |[l6l : [l7ll . 

The paper is organized as follows. In Sec. [2] a generalization of the Skorniakov-Ter-Martirosian equation 
to an arbitrary-angular momentum is presented. In Sec. |3] we show the results of the calculation of the particle- 
dimer scattering lengths for various angular momenta, and the binding energy of the universal trimer states. 
In the Appendix, we give the detailed derivation of the generalized Skorniakov-Ter-Martirosian equation for 
high angular momenta. 



2 Generalized Skorniakov-Ter-Martirosian integral equation for arbitrary angular momenta 

We consider a system of two identical particles (bosons or fermions) strongly interacting with a third particle 
of a different kind. We assume that this interaction is nearly resonant, and thus characterized by an s-wave 
scattering length a, which is much larger than the range of the interaction. On the other hand, we assume that 
the interaction between the identical particles is not resonant, and therefore can be neglected. This is even 
more justified for identical fermions, which do not interact in the s-wave channel. As a result, there is only 
one length scale in the problem, given by a, which we take as the unit length: a = 1 . 

This three-body problem can be solved by an integral-equation approach. The s-wave scattering length 
between the particle and the dimer, in which one of the particles is identical with the unbound one, can be 
obtained by solving the integral equation derived by Skomiakov and Ter-Martirosian QJJ] , and found to be 
1.18 in the equal-mass case. Recently, the calculation was also extended to the unequal-mass case [Q]. Here, 
we derive a similar integral equation to calculate the scattering lengths and the binding energy for arbitrary 
angular momenta. 

The quantity characterizing the particle-dimer scattering in the Z-th wave angular momentum channel is 
the Z-th wave scattering lengths a/, which are defined as a low-energy contribution of the Z-th wave phase shift 
8, 

k 2,+l cotd, = -- + 0(k 2 ). (1) 

Since the phase shift 5; is related with the Z-th wave angular momentum component of the particle-dimer 
scattering amplitude // (k) as 

= I X V (2) 
k cot 0/ — ik 
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Table 1 Mass ratio at which the universal bound state starts to appear, the width of the particle-dimer resonance, and the Efimov 
threshold are shown for various angular momenta and statistics. For the case of / = 1, for example, the universal bound states 
appear at the mass ratios of 8. 172 and 12.917, which reproduces the results of Ref. 1 14]. For higher-angular momenta, we refine 
the values of Ref. 1191 . We also present the Efimov thresholds, above which the non-universality appears and the Efimov states 
can be formed (6112011. 



/ 


statistics 


universal bound state 


width of the resonance 


Efimov threshold 


1 


fermions 


8.172 
12.917 


9.2 
8.1 


13.607 


2 


bosons 


22.637 
31.523 
37.766 


5.9 
11.9 
8.1 


38.630 


3 


fermions 


43.395 
56.166 
67.336 


1.3 
3.7 
5.4 


75.994 


4 


bosons 


70.457 
87.027 
115.534 
102.486 
124.167 


0.14 
0.50 
0.99 
1.27 
0.81 


125.765 



a/ represents the low-energy scattering property [18]. The partial wave scattering length ai can be calculated 
from the following integral equation, whose derivation is given in the Appendix: 




(3) 



aM), 



where ± refers to bosons (+) and fermions (-), a is the mass ratio between one of the two identical particles 
and the third particle, £ is the total energy in units of the dimer's binding energy, Q\ denotes the Legendre 
function of the second kind, and ai (q) is the l-th wave momentum-dependent scattering length whose value at 
q — corresponds to the Z-th wave scattering length a/. By solving this integral equation at the dimer-particle 
continuum threshold £ = — 1, we obtain the /-th wave scattering length. On the other hand, we can find the 
binding energy of the trimer states by searching for the value of £ at which aj diverges. This can be done 
by finding the value of £ at which the eigenvalue of the right-hand side (seen as an operator acting on ci[(q)) 
vanishes. 

We note that in the numerical calculation, we take a finite momentum cutoff instead of the infinite upper 
boundary of the integral in Eq. For a mass ratio smaller than the Efimov threshold listed in Table [TJ we 
obtain convergent results if we take a sufficiently large momentum cutoff. However, if the mass ratio gets 
closer to the Efimov threshold, the cutoff required for the convergence to the universal results becomes larger 
and larger. Above the Efimov threshold, the results depend on the cutoff, as it plays the role of the short-range 
parameter needed to characterize Efimov physics Ill9h . 

We display the results for both bosons and fermions in Fig. Q]to Fig. [5] For the calculations performed 
at a mass ratio smaller than the Efimov threshold given in Table [TJ the results should not depend on the 
momentum cutoff in the integration if we take it as sufficiently large. This is illustrated in the figures by a set 
of curves corresponding to different momentum cutoffs. Below the Efimov threshold, these curves essentially 
superimpose and correspond to the universal curve. Above the Efimov threshold, there is no universal curve 
and the results explicitly depend on the cutoff in this region. 



3 Results of numerical calculations 

In the case of s-wave and bosons (Fig. [TJ (a) and (c)), we have an infinite series of resonances due to the 
formation of bound states. Because there is no centrifugal repulsion and the Efimov attraction prevails for all 
mass ratios, all these bound states are Efimov states, and therefore depend on the cutoff of the integral. As we 
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Fig. 1 S-wave scattering length and binding energy for the bosonic and fermionic cases. (a),(b) .s-wave scattering length, and 
(c), (d) binding energy of trimers are shown. The particle-dimer scattering length is displayed in units of the particle-particle 
scattering length a, while the energy is shown in units of the dimer binding energy. The gray and white regions denote the 
universal and non-universal regions, respectively. The dashed curves correspond to a momentum cutoff of k c = 2000 and the 
solid curve k c = 4000. 



change the cutoff, this series of resonances and bound states gradually shift, while preserving their qualitative 
structure imposed by discrete scale invariance. 

For the .s-wave fermionic case (Fig.Q](b) and (d)), on the other hand, Pauli's exclusion principle suppresses 
Efimov's attraction, and the results are universal and do not depend on the cutoff. The scattering length curve 
presented here reproduces the results already presented in the previous studies fit ul: Bill. 

For the p-wave fermionic case (Fig.[2](b) and (d)), we can see two universal p-wave resonances followed 
by other non-universal p-wave resonances. The two universal resonances are due to the formation of the bound 
states originally found by Kartavtsev and Malykh [14]. The non-universal resonances occur at mass ratios 
larger than the Efimov critical value 13.6, and are thus due to a series of Efimov states. For /7-wave bosonic 
case (Fig. [3), on the other hand, the effective potential is repulsive due to the symmetry of the identical 
particles, and thus the curve is universal and there is no bound state. The widths of the resonances in the 
fermionic case are broad, suggesting that the universal resonances may be observed even when the mass ratio 
is not finely tuned. 

In Fig.[3]to Fig. [51 we present the scattering lengths and binding energy in the higher angular-momentum 
channels. We can clearly see that there are several universal bound states for bosons in the even angular- 
momentum channels, and for fermions in the odd angular-momentum channels. This reconfirms the conclu- 
sion of the previous works lfl6l: [l7ll . The mass ratio at which the resonances occur also agrees with the value 
where the universal trimers were predicted to appear lll6l EPTIl . By looking at these binding energy curves, it 
becomes apparent that they can be grouped in two sets. One corresponds to the Efimov series, where states 
are related to each other by a scale transformation. The other one corresponds to the universal trimer states, 
which were shown in Ref. Ill6ll to be related to a universal function of only two parameters. We note that the 
widths of the resonances become narrower in the higher angular- momentum channels. 
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Fig. 2 P-wave scattering volume and binding energy. Notations are the same as those in Fig. [I] 



Bosons Fermions 




10 20 30 40 50 10 20 30 




no bound state 



10 20 30 40 50 



Mass ratio (X Mass ratio OC 



Fig. 3 fl-wave scattering length and binding energy for the bosonic and fermionic cases. Notations are the same as those in Fig. 
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Fig. 4 F-wave scattering length and binding energy for the bosonic and fermionic cases. Notations are the same as those in Fig. 
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Fig. 5 G-wave scattering length and binding energy for the bosonic and fermionic cases. Notations are the same as those in Fig. 
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Table 2 Mass ratios for some of the atomic hetero-nuclear systems 



species 


statistics 


mass ratio 


7 Li-« 7 Sr- 8/ Sr 


fermions 


12.4 


6 Li- 87 Sr- 87 Sr 


fermions 


14.4 


7 Li- 133 Cs- 133 Cs 


bosons 


18.9 


6 Li- 133 Cs- 133 Cs 


bosons 


22.1 


7 Li- l74 Yb- 174 Yb 


bosons 


24.8 


6 Li- 174 Yb- 174 Yb 


bosons 


28.9 



4 Conclusion 

We have studied the particle-dimer scattering for the unequal-mass case for both bosons and fermions. For 
this purpose, we have derived an integral equation for arbitrary angular-momentum channels, which is a 
generalization of the Skorniakov-Ter-Martirosian equation. With this integral equation, we have calculated 
the particle-dimer scattering lengths for various angular-momentum channels, and found resonances due to 
the formation of universal trimers and Efimov trimers. 

Ultracold atomic gases are currently the most likely candidates for the experimental possibility of ob- 
serving the universal resonances and bound states. In these systems, a very large s-wave scattering length 
can be realized by using a Feshbach resonance, providing access to the universal regime where the scattering 
length is the single parameter characterizing the system [22]. As for controlling the mass ratio, it is possible 
to prepare hetero-nuclear systems with mass imbalance [11], and thus select the appropriate mass ratio from 
a broad variety of atomic species. Some possible atomic combinations are listed in Table [2] We also note that 
several works have also shown that reduced dimensionality or strong confinement favours the appearance of 
similar resonances and trimer states at lower mass ratios lll5l : l23l : l24ll 

If these resonances can be achieved experimentally with ultracold gases, they should provide better control 
over the scattering property than the Efimov resonances. Because atomic Efimov states are inherently unstable 
due to a recombination to lower vibrational states, the atom-dimer Efimov resonances are dominated by 
the strong inelastic collisions and smoothed out. On the other hand, the centrifugal barrier of the universal 
trimers prevents them from undergoing such recombination processes, and makes them stable 11511 . Thus, 
they are promising for controlling interactions in ultracold gases. In particular, the comparatively broad p- 
wave resonances seem to be most feasible. 



Acknowledgements We thank D. Petrov for introducing us to the work of Ref. 11411 ■ We also thank N. Sakumichi and Y. 
Tanizaki for fruitful discussions. 



A Derivation of the higher-angular-momentum Skorniakov-Ter-Martirosian equation 

We here derive an equation to determine the higher angular-momentum scattering length denoted as a/. We start from the three- 
body Dyson's equation 1 25] 



T 3 ( Pl ,p 2 \P) = ±G t (P-pi-p 2 )±iJ ^0G,(P-p i -q)G i (q)T 2 (P-q)T 3 (q,p 2 \P) 



(4) 



where T 2 (q) and T 3 (pi,p 2 \P) are the 2-body and 3-body T-matrices in the vacuum, and Gj(q) is the single particle Green's 
functions in the vacuum (s = 1,2) for identical particles with mass m,, and G,(q) for the third particle with mass m,. The plus 
sign is for bosons, and the minus sign for fermions. Note that all momentum indices are written in four-momentum notation. 

We focus on the on-shell process for p 2 and collision at the total energy P = Pg = {^-,0}, where pi = ™'™ mi is the reduced 
mass between the different particles, and e is the total energy in units of the dimer's binding energy. Then, Eq. © reduces to 

73(pf ,Pa*|flO = ±G t (P -p? -75?) ± J J^G^P-pT 1 -q i )T 2 (Po-q i )T 3 (q i ,p 2 4 \P ), (5) 

-1 

where the overline stands for the on-shell value for the one of the identical particles: q' = {■ 2 ^ m -,q}- Our main motivation is to 
expand this equation into partial waves. To this end, we will first focus on the terms such as 

-G,{P-pi i -pi , ) = — — — - = — — , (6) 

£ I Pi , pj , [pi+PiY e Pj | Pi I PiPi 

2jl ~ 2m,- t 2m; ~ 2m, 2)1 ~ 2jl ~ 2jl ~ m, 
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Fig. 6 Diagrammatic representation of three-body Dyson's equation l(4j. The single curves and the double curves denote single- 
particle Green's functions and two-body T-matrices, respectively. 



which couples two ingoing and outgoing momentum. If this term can be expanded as 

- G, (P - p? - Pi) = £ gi (pi , Pi)P, (cos 6i 2 ), (7) 
z=o 

where 612 is the relative angle between p\ and p2, then we can easily expand Eq. <(5}. In fact, if we expand the T-matrix as 

T 3 tf,p?\Ih) = '£t,(puP2)Pi{cose 12 ) (8) 
;=o 

and substitute this into Eq. (5j, we obtain 

h (j>l , Pi) = Tgl (Pi ,Pz) T (2t + \) f 1 ld <iSi{p\,ii)T2{Po-^)tM,P2), (9) 
where we have used the following formula: 



An 



(10) 



Thus, we only have to solve the above integral equation to obtain the l-th partial wave scattering length. For this purpose we 
have to calculate gi(p\,P2)- From the orfhonormal relation of Legendre polynomials 



J dxP„(x)P m (x) 



2n+l 



we can calculate gi(p\,P2) from the following integration 

2/+1 n 



gl 



(Pi ,P2) = J i dxPi (x) - 



' 2n ~ r 2n ' 2n~ m, 



Here we use the following relation, which we will show below: 
where Q„ (x) is a Legendre function of the second kind. Then, 



gi(pi,P2) = (2i+i)(-iy—Q, 

P\P2 



P\P2 



(ID 



(12) 



(13) 



(14) 



We are only interested in the low-energy value of ti{pi,P2) under p\ = P2, since we are interested in the low-energy elastic 
collision. In fact, the r-matrix and the scattering amplitude are related as 



w,pi-'n) = -— f(pi,P2) 

niDA 



(15) 



where moA = ""'^m.+ 'm^ ' s me reduced mass between the particle and the dimer. This scattering amplitude and the /-th-wave 
scattering length are related as 



fiPUPl) = £(2/+l)Ji(pi)fl(c08fli 2 ) (|pi| = \P2\ 



fi{k) = — 



k 21 

-nk 2 - 



(16) 
(17) 
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where ij is a spherical modified Bessel function and Qi is the Legendre function of the second kind. We have used the following 
formula for these functions in the second line and the third line: 

J 1 dxP n (x)esp(-bx) = 2(-\fi n {b) (b > 0), (32) 

f 'dti H (t)e-» = Q„(p) (p>0). (33) 
Jo 

□ 
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